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Abstract. We present results from a virial analysis of fully nonlinear two- 
dimensional (2D) simulations of the ISM. We discuss the Eulerian Virial Theorem 
in 2D, and describe preliminary results on the virial budget of clouds in the 
simulations. The clouds are far from a static equilibrium, and the Virial Theorem 
is dominated by the time-derivative terms, indicating the importance of flux 
through the cloud boundaries and mass redistributions. A trend towards greater 
importance of the gravitational term at larger scales is observed, although a few 
small clouds are strongly self-gravitating. The magnetic and kinetic terms scale 
f*^ linearly with each other. 
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Or I INTRODUCTION 
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Vazquez-Semadeni et al. (1995, hereafter Paper I) and Passot et al. (1995, 
hereafter Paper II) have produced a numerical model of the interstellar 
medium (ISM) including enough physical agents as to render it feasible to 
perform statistical studies of the clouds formed in the simulations. The simu- 
lations include self-gravity, magnetic fields, parameterized cooling and diffuse 
heating, the Coriolis force, large-scale shear, and localized stellar energy input. 
In the present work, we discuss the Virial Theorem (VT) as it applies to the 
simulations, and present preliminary statistical results from a two-dimensional 
(2D) simulation with a resolution of 800 x 800 grid points, performed specif- 
ically for this analysis. In § II we discuss the VT, applying the formalism 
developed by McKee & Zweibel (1992) to the 2D case. In § III we describe 
the cloud-identifying algorithm and show preliminary statistical results, and 
in § IV we present some remarks and future work. 
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II VIRIAL THEOREM IN 2D 

The VT is obtained by dotting the momentum equation (eq. (lb) in Paper I) 
with the position vector x and integrating over volume. McKee & Zweibel 
(1992) have discussed an Eulerian form of the VT, which is most appropriate 
for our simulations, since they are performed with an Eulerian code. Because 
the simulations are 2D (in order to reach a sufficiently large resolution), we 
must consider the VT in 2D as well. It reads: 

1 d 2 I ( \ 1 rf$ 
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where Tkin = 1/2(J pu 2 dV — § x XipUiUjfijdS) is the kinetic term, T- mt = J PdV— 
1/2 § s PxihidS is the thermal term, M = l/8n § s XiTijhjdS is the magnetic 
term, W = J %ip gidV is the gravitational term, E COT = 2 / Xj(f2 x u)idV is the 
Coriolis term, $ = § s pUir 2 hidS is the flux of moment of inertia through the 
surface S, and p, u, P and ^ are the density, velocity, thermal pressure and 
self-gravitational acceleration, respectively. Because of two- dimensionality, we 
must replace volumes by areas and surfaces by contours in (1). However, we 
retain the above notation for generality. Since in 2D V • x = 2, in equation 
(1) we note the three following interesting points: a) Although magnetic fields 
are present in the surface term M = j s Xi hj dS, (where the Maxwell stress 
tensor is defined as = l/4ir[BiBj — l/2B 2 Sij}), the "classical" magnetic en- 
ergy term E mSLg = l/8n J B 2 dV does not enter the virial equation, so it does 
not provide support against gravity, b) The internal energy E int = / PdV 
does not contain the 3/2 factor as in 3D. Nevertheless, in 2D this term still 
coincides with the total internal energy, because there are only two transla- 
tional degrees of freedom, c) Additionally, it can be shown that the grav- 
itational term / Xi p gi dV does not coincide with the gravitational energy 
-E'grav = 1/2 /p dV as it does in 3D for isolated clouds. Essentially, this is 
due to the slower distance dependence of the gravitational potential in 2D. 

Ill PRELIMINARY STATISTICS 

In order to calculate the terms in equation (1), we have performed a 2D 
simulation similar to the one called "Run 28" in Paper II, but with a res- 
olution of 800 x 800 grid points. In this run we analyze the data shortly 
after turning off star formation, in order to allow for the largest possible den- 
sity gradients (see Vazquez-Semadeni, Ballesteros-Paredes & Rodriguez 1997, 
hereafter Paper III) while still retaining the structure induced by the stellar 
energy injection. We have developed a numerical algorithm to identify clouds 
and evaluate within them the various terms entering the VT, as well as their 
velocity dispersion and mean density. We define a cloud as a connected set of 
pixels whose densities are larger than an arbitrary threshold p t . Previous cal- 
culations (Paper III) have shown that the simulations exhibit similar scaling 
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FIGURE 1. In all panels, the solid line is the identity, (a) l/2\\d$/dt\\ vs. l/2\\d 2 I/dt 2 \\. 
Their near equality shows that the term l/2d<f>/di dominates the virial sum, indicating the 
importance of the variability of the mass flux through the clouds' borders for the total virial 
balance, (b) Volume-plus-surface kinetic terms vs. the virial sum neglecting the l/2d$/dt 
term. The near equality of both terms indicates the dominance of the kinetic terms over the 
remaining ones. This effect may be due to cloud bulk motion and should be eliminated by 
using an instantaneously-at-rest frame of reference for each cloud, (c) Volume vs. surface 
terms for internal energy (pressure) and (d) kinetic energy. The surface terms are seen to 
be comparable to the volume terms in general. The few points with large scatter in (c) are 
likely to correspond to regions of anomalous pressures due to recent star formation, (e) The 
gravitational term W vs. the sum of the remaining virial terms. A trend towards greater 
importance at larger scales is seen. However, a few points at near balance with gravity are 
seen at all scales, (f) Magnetic term M vs. the sum of the kinetic terms. An almost linear 
relation is observed. This is consistent with equipartition between kinetic and magnetic 
modes, if an offset is present, again due to the fact that clouds may have bulk velocities 
with respect to the integration volume. 
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properties as those observed in real interstellar clouds (Larson 1981), except 
for the density-size scaling relation, supporting the possibility that it may be 
the result of an observational effect (see also Larson 1981, Kegel 1989, Scalo 
1990). With this motivation, we have now performed evaluations of the vari- 
ous terms in the VT. We have the following preliminary results: 1.- Both the 
second derivative of the moment of inertia and the last term in the equation 
(1) are dominant in the overall virial balance (fig. la). 2.-Comparing the 
remaining terms, the turbulent terms are seen to dominate (fig. lb). 3. -The 
surface terms (which are often neglected under the assumption of vanishing 
fields outside the clouds) are in general of magnitude comparable to that of 
the volumetric ones (figs, lc and d). 4. -The gravitational term is most im- 
portant at large scales (fig. le). However, there are a few small (low energy 
content) clouds which have large values of the gravitational term. These may 
be the best candidates for collapse and star formation. Their scarcity appears 
consistent with the low efficiency of star formation. 5. -The magnetic term 
and the sum of the kinetic terms are proportional to each other (fig. If). This 
suggests there is equipartition between kinetic and magnetic modes, except 
for a constant factor, which may be due to the fact that clouds have bulk 
velocities with respect to the integration volume. 

IV FINAL REMARKS 

The dominance of the time-derivative and kinetic terms indicates the impor- 
tance of flow through the volume boundaries, contrary to the cases considered 
by McKee & Zweibel (1992). In order to minimize this effect, it appears nec- 
essary to consider Eulerian volumes instantaneously at rest with respect to 
the center of mass of the clouds. However, preliminary attempts suggest that 
the flow through the boundaries cannot be eliminated completely, since the 
clouds are extremely amorphous and change shape rapidly. This work will be 
reported in a future paper (Ballesteros-Paredes & Vazquez-Semadeni 1997, in 
preparation). 
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